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Double negative differential thermal resistance induced by the nonlinear on-site 
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We study heat conduction through one-dimensional homogeneous lattices in the presence of the 
nonlinear on-site potentials containing the bounded and unbounded parts, and the harmonic in- 
teraction potential. We observe the occurrence of double negative differential thermal resistance 
(NDTR), namely, there exist two regions of temperature difference, where the heat flux decreases 
as the applied temperature difference increases. The nonlinearity of the bounded part contributes 
to NDTR at low temperatures and NDTR at high temperatures is induced by the nonlinearity of 
the unbounded part. The nonlinearity of the on-site potentials is necessary to obtain NDTR for 
the harmonic interaction homogeneous lattices. However, for the anharmonic homogeneous lattices, 
NDTR even occurs in the absence of the on-site potentials, for example the rotator model. 
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Heat conduction in low-dimensional systems has be- 
come the subject of a large number of theoretical and 
experimental studies in recent years [l|. The theoretical 
interest in this field lies in the rapid progress in probing 
and manipulating thermal properties of nanoscale sys- 
tems, which unveils the possibility of designing thermal 
devices with optimized performance at the atomic scale. 
As we all know, devices that control the transport of elec- 
trons, such as the electrical diode and transistor, have 
been extensively studied and led to the widespread ap- 
plications in modern electronics. However, it is far less 
studied for their thermal counterparts as to control the 
transport of phonons (heat flux), possibly by reason that 
controlling phonons is more difficult than controlling elec- 
trons. Recently, it has been revealed by theoretical stud- 
ies in model systems that, such as electrons and photons, 
phonons can also perform interesting function 2], which 
shed light on the possible designs of thermal devices. 
For example, heat conduction in asymmetric nonlinear 
lattices demonstrates rectification phenomenon, namely, 
the heat flux can flow preferably in a certain direction 
OHlIj. Remarkably, a thermal rectifier has been exper- 
imentally realized by using gradual mass-loaded carbon 
and boron nitride nanotubes fll| . The nonlinear systems 
with structural asymmetry can exhibit thermal rectifica- 
tion, which has triggered model designs of various types 
of thermal devices such as thermal transistors [8|, ther- 
mal logic gates |l2j . and thermal memory [13j. It is worth 
pointing out that most of these studies are relevant to 
heat conduction in the nonlinear response regime, where 



the counterintuitive phenomenon of NDTR may be ob- 
served and plays an important role in the operation of 
those devices. 

NDTR refers to the phenomenon where the resulting 
heat flux decreases as the applied temperature difference 
(or gradient) increases. It can be seen that a comprehen- 
sive understanding of the phenomenon of NDTR, which 
is lacking at the moment, would be conducive to further 
developments in the designing and fabrication of ther- 
mal devices. The existing studies on NDTR have been 
on models with structural inhomogeneity, such as the 
two-segment Frenkel-Kontorova model [8|, [l4| , the weakly 
coupled two-segment c/) 4 model [l5[ , and the anharmonic 
graded mass model [f|. However, structural asymmetry is 
not a necessary condition for NDTR. In the nonlinear re- 
sponse regime, NDTR can occur in absolutely symmetric 
structures where there exists nonlinearity in the on-site 
potential of the lattice model However, it is still not 
clear whether multiple NDTR can occur in symmetric 
structures. In this Brief Report, we study the exhibition 
of double NDTR in the absolutely symmetric structures 
and find the occurrence of double NDTR. Furthermore, 
we also find that NDTR can also occur in the coupled 
rotator model where the on-site potential is absent. 

In this study, the homogeneous lattice models are each 
described by a Hamiltonian of the form 



H = 



N 

Y — 

^ 2m 

i=l 



Pi 



+ V(xi - x l+ i) + U(xi), 



(1) 
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where pi is the momentum of the ith particle, and Xi its 
displacement from equilibrium position, m is the mass 
of the particles. V(x) is the nearest-neighbor interaction 
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potential, the harmonic potential is used, 

V(x) = \kx\ (2) 

where k is the coupling constant. As for the on-site po- 
tential U(x), we consider two cases shown in Fig. 1. For 
case A (</> 4 model) [l[ 

U( X) = -« X *+* X \ ( 3 ) 

and for case B 

U{x) = -^cos{2^x) + - A x\ (4) 

where a, Uq, and A are the parameters that control the 
shape of the potential. The on-site potential contains the 
bounded and unbounded parts. 

As to obtain a stationary heat flux, the chain is con- 
nected two heat baths at temperature T + and T_, re- 
spectively. Fixed boundary conditions are taken xq — 
xjv+i = 0. For each of the one-dimensional lattice mod- 
els, the equation of motion takes the form 

dH ■ , 

mxi = -- liXi+t,i, (5) 

OXi 

where 7; = 7(^,1 + S iyN ) and & = £+<Sj,i + The 
noise terms £±(t) denote a Gaussian white noise that has 
a zero mean and a variance of 2jksT±, where 7 is the 
the friction coefficient and ks is Boltzmann's constant. 
The dot stands for the derivative operator with respect 
to time t. The local heat flux is given by ji = (x\F(xi — 
x i-i)), where F = — ^ and the notation (...) denotes a 
steady-state average. The local temperature is defined as 
Ti = (mxi 2 ). After the system reaches a stationary state, 
ji is independent of site position i, so that the flux can 
be denoted as j. In our simulations, the equations (5) of 
motion are integrated by using a second order stochastic 
Runge-Kutta algorithm[17] with a small time step (h = 
0.001). The simulations are performed long enough to 
allow the system to reach a nonequilibrium steady state 
in which the local heat flux is a constant along the chain. 

Figure 2 (a) shows the dependence of the heat flux j 
on temperature difference AT for <fi 4 model with a > 
(case A). It is found that there exist two regions of AT, 
in which the larger temperature difference the less heat 
flux through the system, namely double NDTR occurs. 
The presence of a nonlinear on-site potential facilitates 
the occurrence of phonon-lattice scattering, which gener- 
ally becomes more significant for increasing temperature 
and can therefore contribute to a decrease in the ther- 
mal conductivity. For low temperatures (small tempera- 
ture difference), the bounded part of the on-site potential 
dominates the system. For this case, the phonon-lattice 
scattering is important only at sufficiently low tempera- 
tures where the dynamics of the particles is much influ- 
enced by the bounded on-site potential. As the applied 
temperature difference AT increases from zero with T_ 




-200' — 1 1 1 1 ' 

~10 -5 5 10 

X 




-5' 1 1 1 1 1 1 1 1 

-4-3-2-101234 

X 

FIG. 1: (Color online) The shape of on-site potential, (a) 

Case A: U(x) = -f x 2 + ±x 4 (a > 0), (b) case B: U(x) = 
jMg cos(2^) + \x 4 . 



being fixed, the increase in the thermodynamic driving 
force will drive an increase in the heat flux. At higher 
values of AT, however, the effect of phonon-lattice scat- 
tering becomes so significant that the first NDTR occurs. 
But with a further increase in AT, the average tempera- 
ture of the system has become sufficiently high such that 
the particles can overcome the bounded part of on-site 
potential, the phonon-lattice scattering becomes not sig- 
nificant and the first NDTR disappears. At the same 
time, the dynamics of the particles is dominated by the 
unbounded part of the on-site potential. As increasing 
the temperature, the increase in phonon-lattice scatter- 
ing is reflected by the power-law decrease of the thermal 
conductivity and the second NDTR occurs. Therefore, 
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FIG. 2: (Color online) Heat flux j as a function of temperature 
difference AT for case A and B. (a) For case A: a — 10, 
A = 2.0, and T_ = 0.01; (b) for case B: U = 10.0, A = 0.2, 
T_ = 0.001, and N = 32. The insets in (b) give the enlarged 
views of the NDTR. The other parameters are k = 1.0, N = 
32, T+ =T- + AT. 



we can obverse the double NDTR as the temperature 
difference increases from zero with T_ being fixed. Form 
Fig. 2 (b) , we can see that the similar double NDTR can 
also occur for case B. 

Since the results for case A are similar to that for case 
B, we will only focus on case B for investigating the pa- 
rameters dependence of double NDTR. The on-site po- 
tential of case B is composed of two parts: the bounded 
part 



- J^U cos(27rx) and unbounded part jx 4 



Figure 3 (a) shows the heat flux j as a function of AT 
for different values of Uq with A being fixed. For de- 
creasing Uo, the first NDTR region becomes smaller, and 
finally disappears. Therefore, the bounded part of the 
on-site potential contributes to the occurrence of the first 
NDTR. From Fig. 3 (b) , one can find that the occurence 
of the second NDTR is induced by the unbounded part 
of the on-site potential. Figure 3 (c) shows that the two 
NDTR regimes generally become smaller as the system 
size N increases, and eventually vanish in the thermody- 
namic limit. 

Figure 3 (d) shows the effective thermal conductivity 



(c) 

/ 


qOO° QQ *** 

-•-32 

-*-256 


„oO° °oooo£L' 
„o° ******** *f 

** ? 
••••••• 


- —512 


• 




0.01 0.1 1 


10 100 


AT 






FIG. 3: (Color online) (a) Heat flux j as a function of AT for 
different values of Uo at k = 1.0, A = 0.2, and N = 32; (b) 
Heat flux j as a function of AT for different values of A at 
k = 1.0, Uo = 10.0, and N = 32; (c) The total heat flux Nj 
as a function of AT for TV = 32, 256, and 512 at U = 10.0, 
k = 1.0, and A = 0.2; (d) re(AT), AT, and heat flux j as a 
function of AT at U = 10.0, k = 1.0, A = 0.2, and N = 32. 
The other parameters are T_ = 0.001, T+ = T_ + AT. 
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k(AT), heat flux j, and AT as a function of AT at Uq = 
10.0, k = 1.0, A = 0.2, and N = 32. As we know j = 
k(AT)AT. It is found that rc(AT) oc AT -1 - 17 for the 
first NDTR region, and k(AT) oc AT" 1 - 15 for the second 
NDTR region, resulting in j oc T~ 017 for the first NDTR 
region, and j oc T~ - 15 for the second NDTR region. 
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FIG. 4: (Color online) The total heat flux Nj as a function of 
temperature difference AT. The other parameters are T_ = 
0.1, T+ = T_ + AT. 



Note that there is no definite relation between the 
shape of the on-site potentials and the occurrence of 
NDTR. For example, if the term x 4 in case B is replaced 
by the harmonic term (linear potential) x 2 , the shape of 
the on-site potentials does not change essentially, while 
NDTR at high temperatures will disappear. Therefore, 
it is the nonlinearity, not the shape of the on-site po- 
tentials, that determines the occurrence of NDTR. The 
nonlinearity in the on-site potentials is necessary to ob- 
tain NDTR for the harmonic interaction homogeneous 
systems. Now we return to the anharmonic homogenous 
systems and check if NDTR can occur in the absence of 
the on-site potentials. From the previous work [l6j . one 



can find that NDTR can not occur in pure harmonic and 
Fermi-Pasta-Ulam model. However, in that work, the 
rotator model was not considered. The simplest exam- 
ple of rotator model with nearest neighbor interactions 
lies in the class (1): V(x) = 1 — cos(x) and U(x) = 0. 
From Fig. 4, we can see that NDTR occurs in the ro- 
tator model, which is related to excitation of nonlinear 
localized rotational modes of the chain [18|. In addition, 
NDTR can also occur in the anharmonic graded mass 
lattices Therefore, for the anharmonic systems, the 
on-site potential is not a necessary condition for NDTR. 

In conclusion, we study heat conduction through the 
onc-dimcnsional homogeneous lattices with the nonlin- 
ear on-site potentials. The on-site potentials are com- 
posed of two parts: the bounded and unbounded parts. 
From nonequilibrium molecular dynamics simulations, 
it is found that double NDTR occurs as the tempera- 
ture difference increases. The occurrence of NDTR at 
low temperatures is caused by the nonlinearity of the 
bounded part, while the nonlinearity of the unbounded 
part induces the occurrence of NDTR at high tempera- 
tures. In addition, we also find that NDTR even occurs 
in anharmonic homogenous lattices without on-site po- 
tentials, for example the rotator model. Therefore, the 
on-site potential is not necessary to obtain NDTR for 
the anharmonic lattices. It is also found that the regime 
of NDTR becomes smaller as the system size increases, 
and eventually vanishes in the thermodynamic limit. The 
observation of double NDTR in homogeneous systems 
shows that the nonlinearity of the on-site potentials is 
very important for designing NDTR devices. It is possi- 
ble to design the thermal devices with the more complex 
functions by using the occurrence of double NDTR. 
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